Persistent current in quantum torus knots 
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I consider the quantum interference of electrons moving along knotted trajectories under external 
magnetic field. The induced persistent current is formulated in terms of characteristic parameters 
that classify the torus knot geometry. The current is found to show a periodic oscillation whose 
period depends strongly both on the knot's class and the field direction. The shift in the oscillation 
period caused by geometric distortion is also discussed. 



State-of-the-art fabrication techniques have enabled 
to realize microscopic "knotted" objects based on or- 
ganic and/or inorganic materials. A knot means a self- 
entangled closed loop that cannot unraveled except by 
cutting the loop; Figure [T] illustrates the simplest exam- 
ples of knots. The presence of knotted molecules was 
first identified in DNA p], [1[ , and then naturally occur- 
ring proteins 0, 0] • After the field of chemical topology 
blossomed, diverse molecular knots endowed with differ- 
ent types of topology have become in reality [5j. Added 
to the molecular-scale entities, micrometer-scale knots in 
carbon nanotube bundles [6] and those made of triblock 
copolymer Q have also been realized. To elucidate the 
knotted effects on their physical properties, the isolation 
and precise characterization of their shapes are required. 
However, the task is still a challenge even today, which 
is why little has been reported for the shape-property 
relation in knotted materials 

An important consequence of the knotted structures, 
provided they are metallic, is thought to arise in magnetic 
response of internal mobile carriers. Generally when a 
closed loop of a quantum wire is threaded by a mag- 
netic flux <£, a persistent current will be elicited due to 
quantum interference of electrons along the loop. For 
instance, a micrometer-diameter ring can support a per- 
sistent current of / ~ 1 nA at temperatures T < 1 K 
and the magnitude of the induced current oscil- 
lates with increasing with a period of the flux quantum 
<l> = h/e. This scenario may hold true no matter if the 
loop is knotted or unknotted. For knotted cases, however, 
the flux is tangled in self-intersecting, multiple surfaces 
surrounded by the self-avoiding closed loop; hence, the 
resulting persistent current possibly depends on the field 
direction and the knotted geometry. Evaluating these 
dependences will facilitate the development of molecular 
machines and other bottom-up nanodevices that operate 
under magnetic field. 

In this Brief Report, I describe the persistent current 
phenomena arising in a certain group of knots, called the 
torus knots. Persistent currents that occur in general 



torus knots are formulated as a function of geometric pa- 
rameters of the knots, which reveals the strong anisotropy 
of the current magnitude with respect to the direction 
of the external magnetic field. Analytic expression for 
the period of the current oscillation with increasing the 
field strength is also derived, showing the way how the 
knotted geometry affects the oscillation period. These 
findings will provide the first step toward quantitative 
analyses on the magnetic response of realistic molecular 
knots and other micrometer-scale knotted conductors. 

To proceed analytic discussion, we employ a one- 
dimensional coherent electron system having a torus knot 
geometry. Torus knots consist of a special group of knots, 
satisfying the condition that they all lie on a donut- 
shaped surface. The class of a torus knot is uniquely 
identified by a pair of relatively prime integers p and q 
[l2j ; the one integer p specifies the number of times the 
curve wraps around the rotational axis of the torus, and 
the other q specifies the number of times the curve passes 
through the hole of the torus. If p and q are not co-prime, 
then we have a collection of two or more identical knots 
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FIG. 1. (color online) Diagram of the four simplest torus 
knots labeled by (p, q) = (2, 3), (2, 5), (3, 2), (5, 2), respec- 
tively. The lower right panel illustrates how the (2,3) torus 
knot wraps around the underlying torus colored in blue. 
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that all lie on the same torus. The lower panel of Fig. 1 
shows schematically the way how a (2,3) torus knot wraps 
around an underlying donut-shaped surface. 

In terms of Cartesian coordinates, the (p, g)-torus knot 
is parameterized by 

x(0) = R + ecos(q6) cos(p#), 



2/(0)= R + ecos(qO) sin(p0), 
z(0) = £ sin(g0) 



(i) 



with a variable (0 < < 2tt). The two constants R and 
e (0 < e < R) are called by the major and minor radii 
of the torus, respectively; e is the radius of the tube that 
consists of the torus, and R is the radius of the circle 
made of the tubular axis. It follows from Eq. (pQ) that 
the (p, —q) torus knot is the mirror image of the (p, q) 
torus knot and the (— p, —q) torus knot is equivalent to 
the (p, g)-torus knot except for the reversed orientation. 
In the following, we assume that p, q > without loss of 
generality. 

Suppose that the (p, g)-torus knot, designated by T, 
is subjected to a uniform magnetic field B. When an 
electron migrates between two points r a and along T, 
it earns an additional quantum phase Xa^b defined by 



A ■ dr. (e > 0) 



(2) 



Here, A = (B x r)/2 is the vector potential associated 
with B. The total phase shift \T f° r the electron going 
around T reads as 



xt(b) 



-2tt 



MB) 



(3) 



where & T (B) = Bdn is the sum of fluxes that thread 
the multiply self-intersecting p surfaces surrounded by T, 
and n is the unit vector normal to the surface whose total 
area is designated by S T . 

We are ready to derive the persistent current I in a 
^-length torus knot containing N electrons with the effec- 
tive mass m* . The current I n carried by a single electron 
in the nth eigenstate is I n = ev n /£ = ehk n / (ra*£), where 
the wavenumber k n is given by 



k n — 



2tt 

T 



$0 



n = 0,±1,±2,... (4) 



The total current I at zero temperature is obtained by 
summing the contributions from all eigenstates with en- 
ergies below the Fermi level. Remind that / for odd TV, 
denoted by 7 dd? differs from that for even TV, denoted 
by /even, In fact, we can prove that 



(N-l)/2 
lodd = 2 X In 
n=-(N-l)/2 



Mg) 

) — -7 

$0 



(5) 



for -1/2 < $ r /$o < 1/2, and 

N/2 
n=-AT/2+l 



Mg) 



(6) 



for < $ r /<I>o < 1 with the definition of J = 
2N eh / (m* i 2 )] Note that 7 dd and I even have the same 
period such that I a ($ T ) = I a (®r + $o) W = odd or 
even] and I a (0) = 0. Since precise control of N is dif- 
ficult in experiments, we assume an ensemble average 
over many realizations of quantum torus knots to obtain 
(^odd + ^even)/2. As a result, we arrive at 



I = -h 



for0<|^<i, (7) 



and I = at $ r = 0, J($ r ) = J($ r + ($o/2)). It follows 
from Eq. ([7j) that every geometric information (except 
for £) is involved in the flux-sum <l> r . The remained task 
is, therefore, to reveal the dependencies of $> T (B) on the 
direction of the field B and the geometry of the knot. 

To deduce the explicit form of & T (or equivalently xt)> 
we use the relation 



A dV AO ( 

dr = Te de = VdB e - 



dO and B = £? ? -e ? -, 



(8) 



B x ,B y ,B z for i = 1,2,3, 



where m = x, z and ^ - ^,^,^2 ™ . - 
respectively, and e$ is the i^-directed unit vector |13j ; 
the summation convention with respect to i was used in 
Eq. From Eqs. © and ©, we find 



Xr(B) 



— e 
2ft 



2tt 



(xy f - yx) 



+B x (yz'-zy') + B y (zx'-xz') d6, (9) 



where x' = dx/dO so do ?/ and 2/. Referring to Eq. (pQ), 
the integrand of Eq. (j9j) is written in terms of a trigono- 
metric series. For instance, we have 

xy'-yx' = pR 2 + 2pRecos(qO) + ^e 2 l+cos(2<?0) , (10) 

and thus straightforward integration yields 

J (xy' -yx')d6 = 2p(7rR 2 + ^s 2 ) : (11) 

where we took into account that p(^ 0) and q(^ 0) are 
relatively prime integers. 

In a similar manner, it is easy to derive 

yz' - zy' 



Rer 

Tl 

J2 



(p - q) sin[(p + q)0) + (p + q) sin[(p - q)0] } 
|4g sin(p0) +p sin[(p- 2g)0] -p sin[(p+2<?)0] | , (12) 
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FIG. 2. Landscape of spatially modulated minor radius s(6) 
of the torus knot. Three typical patterns are presented, all 
of which fluctuate across the center line e — R/2 with the 
amplitude Se = R/4. 



and 



zx — xz 
Re 



= -y { (p - q) cos[(p + q)0) -(p + q) cos[(p - q)6] } 
-^4q cos(pO) +p cos[(p- 2q)0] -p cos[(p+2q)6] |(.13) 



e 

"~4 



Note that p±q ^ and p±2q ^ 0, since p, q are coprime. 
Hence we have 



/ (yz'-zy , )d0= / (zx 7 - arz 7 ) dfl = 0. 
Jo Jo 

Substituting the results into Eq. ((9]), we obtain 
XT {B) = ^xpB z (ntf + ^e 2 ), 



and 



$ r (B)=pB z (nR 2 + ^s 2 ). 



(14) 

(15) 
(16) 
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FIG. 3. (color online) Probability distribution function P(A) 
of the excess term A defined by Eq. (|17[) . The field direction 
is set as B/\B\ = (0, 0, 1) for (a) and B/\B\ = (^=, ^, ^) 
for (b). 



The formula (fT6|) is the main fruit of this article, de- 
scribing the effects of the knotted geometry and field di- 
rection on the periodic oscillation of /(^r)- The follow- 
ing three peculiarities can be drawn. First, the trans- 
verse components of the field, B x and B yi play no role in 
the quantum interference in the present systems. Only 
the vertical component B z determines both the quantum 
phase shift Xr and the oscillation period of /, although 
the other two components should be reckoned in the con- 
tour integral defined by Eq.(j2j). Vanishing the B x - and 
^-contributions owes to the zero-sum rule expressed by 
Eq. (fl4]h which holds for any choice of (p,q). 

Second, I(& T ) m (p, g) -torus knots is independent of 
depending only on p in the form of / ex p. This fact is 
trivial if e —± 0, namely, the knot is reduced to a simple 
i?-radius circle along which the electron travels p times 
for one duration. But for finite £, (p, g)-torus knots devi- 
ate geometrically from the i?-radius circle and those with 
different values of q (but the same value of p) exhibit 
different entangled structures extending over the three- 
dimensional space, which makes the ^-independence un- 
obvious. 



Third, the minor radius e gives to <l> r an excess con- 
tribution on the order of e 2 . For real knotted materials, 
it may amount to e ~ R/2 more or less; hence the ex- 
cess term ire 2 / 2 shown in Eq. ([T6]) causes <l> r to increase 
about ten percent compared with the case of e = 0, no 
matter what value p and q possess. It should be em- 
phasized that all the three peculiarities listed above are 
implications of the single formula ([16]) . which we have 
analytically derived with neither aid of numerical com- 
putation nor mathematical approximation. 

One may be interested in whether geometric distor- 
tion of the knots causes a change in the formula ([T6]) . It 
is certain that the trajectories of actual knotted materi- 
als deviate from the ideal torus knot geometry because 
of mechanical distortion and thermal disturbance. Such 
the imperfection effects can be evaluated in part by con- 
sidering a distorted torus knot whose minor radius e is 
spatially modulated, as demonstrated below. 

Figure [2] shows typical profiles of the ^-dependent 
e(0) we have considered, wherein every curves fluctuate 
smoothly around the centerline e = R/2 with the ampli- 



tude of 5e = R/4. These modulated e{6) are generated 
by the Fourier filtering method; see Refs. [3, [H| for de- 
tails. The spatial modulation in e leads a breakdown of 
the zero-sum rules ([T4|h thus revives the B x - and B y - 
contributions in the integrand of Eq. ((9]). As a result, 
the formula ([T6]) should be revised as 

* T (B) = pS* x ttR 2 [ 1 + A(B, fc) ] . (17) 

The excess term A takes various values depending on 
the e{0) configuration of individual distorted (p, q) torus 
knots. Hence, we evaluate A for various £r(#)-realizations 
and field directions to obtain the probability with which 
a specific value of A occurs. 

Figure [3] demonstrates the probability distribution 
functions P(A) under the conditions of: B/\B\ — 
(0,0,1) for (a)andB/|B| = (^,^,^)for (b). We set 

e = R/2 and (p, q) = (2, 3) for all the plots. It is found in 
Fig. [3ja) that a sharp peak emerging at 5e <C R becomes 
broadened with increasing fe, together with the slight 
shift in the peak position from A ~ 0.13 to A ~ 0.17. 
Similar behavior is observed in Fig.[3^b), while two peaks 
(instead of one) moving apart from each other are ob- 
served. Attention should be paid for that in the latter 
panel, A for large Se reaches minus one in rare setting, 
which results in <I> r = as seen from Eq. ([TT]) . The same 
situation takes place incidentally when the field direction 
deviates significantly from the vertical one. This fact im- 



plies the preference for distortion-free knots preparation 
in order to observe a well-defined periodic oscillation of 
7(<I> r )in experiments. 

As a closing remark, it is worthy to mention that the 
history of persistent currents dates back to the early days 
of quantum mechanics; nonzero orbital angular momen- 
tum of aromatic molecules were already noticed many 
decades ago [H, [I?}. Nowadays, it is revealed that real- 
istic knotted materials may have twisted internal struc- 
tures such as twisted 7r-electron orbitals in torus-knot- 
shaped organic molecules [18|. The twisted nature is 
known to affect the aromaticity of knotted molecules 
[~L9^ or yield zero-field quantum phase shift [20|, which 
we have omitted in this work. It is also interesting to 
note that quantum confinement into bend hollow cylin- 
ders such as deformed carbon nanotubes 0, [2l[ produce 
an effective electro-static potential field; the geometry- 
induced potential can result in a qualitative change in 
sing le- p article states [To|, HH Hil and collective excitations 
[24 [25| of electrons in the system. Though the present 
work is based on a simplified model, I believe that the 
results obtained provoke further quantitative discussions 
as to knotted material properties taking into account de- 
tailed atomic/molecular configuration and the geometry- 
induced field effects commented above. 
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